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Seminary 

The theory of galactic heavy-ion fragmentation has 
been furthered by incorporating a T- matrix approach in- 
to the description of the three-step process of abrasion, 
ablation, and final-state interactions. The connection 
between this T-matrix and the interaction potential is 
derived. The resulting transition rate is shown to be in- 
dependent of the choice of the initial time. For resonant 
states, the substitution of complex energies for real 
energies is formally justified for up to third-order 
processes. 

The previously developed abrasion-ablation fragmen- 
tation theory is rederived from first principles and is 
shown to result from time-ordering, classical-probability, 
and zero-width resonance approximations. Im- 
provements in the accuracy of the total fragmentation 
cross sections would require an alternative to the latter 
two approximations. Since a more rigorous test of the 
theory would be to compare theoretical and experimental 
differential cross sections, a Lorentz-invariant differen- 
tial abrasion-ablation cross section is derived which ex- 
plicitly includes the previously derived abrasion total 
cross sections. This result requires the use of the time- 
ordering and classical probability assumptions. It is 
demonstrated that spectral and angular distributions 
could be easily obtained from the general Lorentz- 
invariant form. Future success in calculating these 
distributions will require the evaluation of the ablation 
T-matrix, which is the remaining formidable task. 

Introduction 

Research efforts are currently under way to develop 
methods for protecting astronauts from the potentially 
harmful effects of cosmic rays. These effects are par- 
ticularly important for astronauts on long-duration mis- 
sions and for career space workers making repeated 
journeys into space in the era of the Space Station and 
the Space Transportation System, because a biologically 
significant component of cosmic rays consists of 
relativistic nuclei (ref. 1). Theories are presently being 
developed to accurately describe their interactions with 
matter (refs. 2 through 6). In previous work, an optical 
model potential approximation to the exact nucleus- 
nucleus multiple scattering series has been used, within 
the context of eikonal scattering theory (refs. 3 and 4), to 
accurately predict total abrasion cross sections. This 
method has been extended (ref. 6) to calculate isotopic 
and elemental total fragmentation cross sections by 
essentially multiplying the total abrasion cross section by 
the following two factors: (1) a charge dispersion fraction 
giving the probability that z of the abraded nucleons are 
protons, and (2) a compound nucleus decay probability 
which describes the de-excitation (ablation) of an excited 
projectile nucleus prefragment. Final-state interactions 


(ref. 7), which describe the interactions between the 
abraded nucleons and the remaining projectile prefrag- 
ment, arc not yet included in the theory. 

Because of greater flexibility in comparing theory 
with experiment, it is advantageous to develop methods 
for calculating Lorentz-invariant differential cross sec- 
tions. Aside from direct comparisons with experimental 
data, angular and spectral distributions are easily obtained 
for any reference frame (e.g., nucleus-nucleus center of 
mass, nucleon-nucleon center of mass, laboratory, etc.). 
Angular distributions obtained from eikonal theory are 
usually evaluated in the laboratory frame, and transfor- 
mations to another reference frame are laborious (ref. 8). 
In addition, eikonal theory is not readily generalized to 
calculate the Lorentz-invariant differential cross sections 
which are often used to present experimental data. 
Rather than simply multiplying the abrasion cross section 
by approximate factors to obtain ablation effects, it is 
desired to include ablation in a more fundamental and 
exact manner by explicitly calculating ablation matrix 
elements. 

With the preceding considerations in mind, a first- 
prinriples derivation of the transition rate is presented 
herein. The transition rate formally includes abrasion, 
ablation, and final-state interaction matrix elements and 
is easily generalized to produce Lorentz-invariant dif- 
ferential cross sections^ The interaction potentials V and 
transition operators T introduced herein are used to 
describe the three-step process of abrasion, ablation, and 
final-state interaction. This generalization is the main 
conceptual difference between their use in references 2 
and 3, where they were directly connected to an optical 
model description of abrasion only, and their use in the 
present work. 

Evaluation of Amplitudes 

The fundamental equations for the probability 

amplitudes are derived in this section. The full Hamilton- 
ian H consists of an unperturbed piece H 0 and an addi- 
tional interaction V y such that 

H=H 0 +V ( 1 ) 

The eigenstates and eigen-energies of H Q are obtained 
from 

H 0 \n> -ej *> (2) 

The full wave function ¥ is determined from the time- 

dependent Schrodinger equation 

H* = « U (3) 

This wave function is expanded in terms of the complete, 



orthonormal set of unperturbed states | n> generated by 
H a as 

* = E C nW ex P ( - * R t/fi ) ! « > (4) 


Substituting equation (4) into equation (3) yields (refs. 9 
and 10) 

= 4- L c n (,) eX P ("wO V kn (5) 
at m n 

which, for a time-independent interaction, has the 
solution 

<*(') = K + ^ c n (/')exp(^ n f')df' (6) 


where 

* <«* - O'" (7) 

and 

K *n * <*i K|,l> (8) 

Equation (S) is exact, is of fundamental importance, and 
does not require V kn to be small, as might be assumed. 

The weak incompletely coupled approximation 
(WICA) is defined as 

dcdt) | 

~dT * In C| (0 exp {iUk>t) Vk > (9) 

where 

<*>*,= 00 ) 

and 

V ki m<k\V\i> (11) 

with the initial, fully populated state at time t 0 denoted 
by |i>. This approximation (eq. (9)) uncouples c k (t) 
from all other amplitudes except c,(/) and suggests that 
the regeneration or depletion of a particular stite 
depends only upon its being fed from or decaying to one 
other state. As is discussed subsequently, this approxima- 
tion is equivalent to neglecting higher-order terms in the 
interaction matrix elements. 

The strong incompletely coupled approximation 
(SICA) sets 

c*< 0 = (12) 


in equation (3) to give 


dc k (t) i 

-^-*-exp fr k flV u (13) 


This equation is valid only for times which are short com- 
pared with the lifetime of the initial state (i.e., 
c k (t) « c,(0). If the perturbation V is transient and 
small in magnitude, the amplitudes c k (t) will remain small 
for all times. Thus, equation (13) is most commonly used 
for scattering problems with small pecturbing interaction 
potentials. In contrast, the WICA (eq. (9)) is appropriate 
when studying a decaying state over a period of time 
equal to several half-lives. Therefore, the WICA is ap- 
propriate for evaluating resonance states, since their 
complex energy widths result in the usual exponential 
decay. The solutions of equations (9) and (13) are iden- 
tical except for the appearance of complex energies in 
equation (9) where real energies are found in equa- 
tion (13). Thus, one of the major results of the present 
work is the demonstration that real energies can be simply 
replaced by complex energies when studying resonances. 

Real-Energy Transition Amplitude 

In this section, equation (3) is integrated directly to 
yield an “exact” amplitude. A convergence factor, 
a— 0 + , will be inserted to render the integral convergent 
for a nonfinite initial time {t 0 — oo), where at is much 
smaller than unity. Thus, equation (6) becomes 

x c„(t) exp Uui kn r + at') dt' (14) 

and analogously, 

'* i 

* / q (','') ex P (*V,"+ at ") dt" ( 15 ) 

‘o 


In equations (14) and (15), the correct method is to take 
the limit as a-0 + after taking the limit as Also, 

the convergence factor a is not required (Set = 0 before 
performing the integrals) for a finite initial time (e.g., 

'• = °>- 
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Inserting equation (IS) into equation (14) and 
expanding yields and iterative result as follows: 


Equation (17) is the exact final-state amplitude derived 
from equation (14). 


c *(') = hi + ^ y kn jTcxp (L> k ,r + at') dt' 
+ at") df"J 

■ "'1 **P ("it/" + at") dt"| 


exp (u kn t' + a/')dt' 


+ at"')d/ 


x exp(iu kn t' + at’)dt' + . . . 


(16) 


Evaluating these integrals gives 

, , , 1 , «P <**t/ + -exp («*,/„ + «g 

c *<'> = h, + * v k, — 

“ -u>.+ #or 


Real-Energy Perturbation Theory and the SICA 

In this section, a much more compact form for the 
amplitude is derived through the introduction of a 
T-matrix. An ansatz suggested (ref. 10) by the SICA 
(eq. (13)) is written as 

c k (t) = & ki + -rjjT k , f ex P («*>*,/' + aOdt' <*8) 


where the convergence factor is again inserted to render 
the integral convergent for nonfinite initial time 
(t 0 — -oe ), and T ki is a transition matrix element whose 
relationship to v k, will be determined. Evaluating the in- 
tegral in equation (18) and taking the limit as t Q — -ao 
before taking the limit as a^0+ yields 

T 

c k^o - -®) = h , + «P<*V + al) < 19 > 


For a finite value of t 0 , the convergence factor is not re- 
quired and equation (18) becomes 




Tk, 

-n* ki 


exp (h> ki tj {exp( ki ki (t - - 


1} 

( 20 ) 


+ 


j_ v (exp (iu k ,t + 2a/) - exp (to t / 0 + 2 a/„) 

iPn k ” n '\ (-<*>„, + »)(-(*>*, + ») 


Substituting equation (18), and the analogous expression 
for c n (t,t 0 oo), into equation (5) yields the relation- 

ship between T and V as 


ex P (^Vo + al o> |?*P (to»n t + ap- exp t iu> tn t 0 + 2a 
(-u> m + «)(-«*„ + ») 


h. = y k , + E 

n 


y k„T m 


+ ia) 


( 21 ) 


+ tt££ ****** 

n n I 

i exp lu k jt + 3a/) - exp (u kn t 0 + 3a g 
(-«>*, + '3a) (-tai /( + «) (-w m + /2a) 

| exp(to m / 0 + 2a I J [exp (iu kn t + at) - exp(iu kn l 0 + ottj) 

(-“to. + «)(“«« + «)(-"„/+ «) 


This same result (eq. (21)) can be obtained from the exact 
solution (eq. (17)) by taking the limit as t 0 —-» in equa- 
tion (17), substituting equation (19) into the left-hand 
side of equation (17), and solving for ^to- 

identifying the denominator in equation (21) as the 
Green function G„, and using a modified Einstein sum- 
mation convention for triply repeated indices enables 
equation (21) to be written as 

T k,= y k)+ y knG m T n , (22) 


e *P (*»*/'<> + af o> [exp(to k S + 2at)-exp(iu kl t 0 + 2atJ] 
(-“*/ + /2a)(-u> /( + to) ( - u) n/ + to) 


or, in operator form, as 

T= V+ VGT 


(23) 


(17) This is the usual result from scattering theory (ref. 10). 
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Because of its extreme compactness and simplicity, 
equation (19) would be useful for the fundamental 
amplitudes in applications to the heavy-ion fragmenta- 
tion problem. There are, however, some questions con- 
cerning its general applicability: (1) how do the solutions 
depend upon the choice of the initial time t Q ; and (2) how 
does one include a description of resonances, which re- 
quire complex energies, rather than the real energies used 
previously. These questions are addressed in the follow- 
ing sections. 

Transition Rates and the Initial Time 

The choice of a nonfinite ir itial time (t 0 — -oo) is ap- 
propriate in a scattering situation to ensure that possibly 
spurious transient terms are not introduced when the 
scattering region suddenly experiences the full potential 
at time / 0 . It is, however, worthwhile to determine what 
effect the choice of t Q has on the experimentally 
measurable transition rate (or cross section). 

Clearly, the amplitudes given by equations (19) and 
(20) differ for different values of / 0 . The transition rate is 
defined as 


w _ lim 
w ~ dt 


(24) 


where, for discrete final states, the probability is 

m =£ p k it) ( 25 ) 

k 

and for continuous final states, 

PV)=/p(f k )P k Wdf k (26) 

with the probability related to the amplitude as 

p k (f) = k*(0i 2 (27) 

From equation (19), 

| -—II 1 - „y ,p(2»qlr„|! <2*> 

which, upon taking the limit as a— 0 + , becomes 


A + J t \Wo ~ -»>r = f 2 T k p (29) 


2 2x t 


Since the right-hand side of equation (29) does not 
change for the limit as t ^oo f the transition rate becomes 


" = f (30) 

Equation (30) is a fundamental result. In a similar man- 
ner, equation (20) yields 


d 

dt 


k*agi 2 


1 IT . „ 2sin(u k ,t) 
* |7 *' 


(31) 


Equation (31) yields 


<32 > 

which is identical to equation (29). Hence, the transition 
rate is independent of the choice of the initial time. 

Resonances and Complex Energies 

In previous work (refs. 6 and 11), an abrasion- 
ablation model was developed to describe heavy-ion 
fragmentation. A possible source of disagreement be- 
tween theoretical predictions and the available ex- 
perimental data might be uncertainties in the excitation 
energies of the excited projectile prefragments (refs. 6 
and 12). In reference 6 it was assumed that the prefrag- 
ment could be treated as an excited compound nucleus. 
An alternative approach would be to use a nuclear 
cascade description (lef. 13). Since a correct description 
of the excited prefragment is essential in predicting 
fragmentation cross section , a description incor- 
porating an intranuclear ca^ade followed by a com- 
pound nucleus evaporation process may be required. 

The treatment of resonances, such as for the com- 
pound nucleus (ref. 14), requires the use of complex 
energies (ref. 9). In first-order theory, Merzbacher 
(ref. 10) has shown how to use complex energies to 
describe resonances, and Norbury and Deutchman 
(ref. 9) have shown the same thing in second-order 
theory. Both of these presentations, however, depend 
upon the use of the WICA (eq. (9)). As previously men- 
tioned, the consideration of final-state interactions, 
together with abrasion-ablation, requires, in principle, 
third-order matrix elements. To accomplish this it is 
shown subsequently how to include complex energies 
beyond second order without the use of an incompletely 
coupled approximation. 

The basic second-order process involves a transition 
from an initial state |/> to an intermediate resonating 
state | /?>, followed by subsequent decay to the final 
state | Ar> . This transition involves the use of three cou- 
pled equations of the type given in equation (5), whereas 
Merzbacher (ref. 10) uses only two coupled equations to 
extract the first-order process. The rate of change of the 
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intermediate-state coefficients is given by 

777 = -jrL C/(Oexp(w rt/ 0 V nl 
at m i 


dcM) V mM \ „ 

- -~r C " (,) + ~ Z C/W exp (*V> (39) 

(33) 

which, upon substitution of equation (37) becomes 


where the rates of change of the c t (t) are given by 

^7^ "~rT, c m ( '> e *P (*"/*'> V lm < 34 > 

dt m m 

In reference 9, Norbury and Deutchman solved these 
equations by assuming the WICA form for equation (34) 
as 

— = - c m (Oexp (iu> lm t) V lm (35) 

dt iff 


*„<0 v. 


dt m 


~ exp (u»0 + 


ill 


HE E I v »/l 2 **P 
W /*» 



/*/! m*n 


x /' 0 m (O e..F (*» lm t’)df 


(40) 


In the present work, the full equation (34) is considered. 
A key factor, however, in the solubility of the three 
coupled equations (eqs. (5), (33), and (34)) is the separa- 
tion of equation (34) as 

~ = jjj c n(0 “P (fa>,„/) y, n 

+ i C *nW eXp V tm < 36 > 


Because of the substitution of the Kronecker delta term 
from equation (37) into the summation term in equa- 
tion (39), n * i in V nr Because /* n in equation (39) the 6 fl 
in equation (37) implies that w*i in l^ w . 

The solution of equation (40) c n {t) appears com- 
plicated but can be simplified by reducing equation (40) 
to an algebraic expression through the use of a Fourier 
transform pair, defined as 


/ oo 

J n { a?) exp ( - locf) r/w (4 1 ) 

equation tJ5), wmen is the wica contribution. 00 


Integrating equation (36) yields 


and 


C/( r) 6 „ + Js f Cn(nexp {L , ln r )d r 

+ T L V ^S‘ ^')«P dt' (37) 


** m*n 'o 


where 


<v(' 0 ) = 


(38) 


Equation (33) is rewritten ai 


/„(«■-) = J* t exp (io;f) dt (42) 

Assuming 

<■„(' £ '„) = 0 (43) 

for ihe intermediate state, equation (42) becomes 


/„(*>) = - f exp <fc 
2 x »/„ 


(iu>f) dt 


(44) 


Multiplying equation (40) by exp [i(«) + ia)t] and in- 
tegrating from t 0 to oo yields 
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dc n (t) 

exp [<(oj + «*)/) dt 

dt 


- h y |v I 2 ^ 

UW /*„ fa n i + o))-a 

V f» 

— I exp (t(w_. + w + /a)/] dt 
at J, » 


J y K«/ ^/rn 

" i(<rf + w w/ )-a 
1 } M/i //M/i nt 

V e* 

+ — J. c„(t) exp (f(w + ia)t] dt 
it! '<> 


* f exp [/(«„,„ + w + ia)f ' ) c m (/‘ ) <//' ( 46 ) 

*o 



Solving for /„(w) and taking the inverse Fourier 
transform yields 

+ *a)/J 

2*7 J f n {ui) exp ( - Kit) du 

- QO 

x fjf' c„(f) exp d/J dt 


= ~ exp [*«„,+ ia)g 

*.i E I 

(,ff) (*n m±n 


exp - f)l du 

x J 

- <*> (w + u> m + ia) l~{o) + /a) + £ (w)] 

fOO 

x J' exp [/(«„, + w + /a)/] 


TfZ Z V nl v lm 

n M/i m±n 

* ex P t^tm 1 )^j dt 

( 45 ) 

/**{/ exp + + C m( / ) <rt'je\pHu/)db 



J_ lfa„i + iai) - <*)(-(« + id) + £(u)l 


( 47 ) 


where the a factors have been inserted to render the in- with 
tegrals convergent. Upon evaluation of these integrals, 
equation (45) becomes 


, . . nrt 
£(<*») ~ — + 

fi 




z 




t*n w m 


U) n / + W + kx 


(48) 


-2 *-/,V + 

V ni exp [Hu ni + u + ‘a)t 0 ) 
fi + w + /a 




Evaluation of the integrals in equation (47) is simplified 
by assuming (ref. 10) that 


£(w) = {(0) 


and inserting Dirac's identity (ref. 15) 

At. i r 

lim £(0) = — - - - 

ot-0 + it it 2 


(49) 


(50) 
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where the energy shift Ae„ and the width T are defined as 

4 *- k - + E ,7^ <51) 

ten * 1 

and 


For f<f 0 , a contour taken in the upper half plane of 
figure 1 gives 

I 0 (t < g = 0 (58) 

since no poles are enclosed. For J Q (t > the lower 
contour, which encloses the poles, is used to yield 


Z I y J «<*«-*/) (52) 

2 ten 

The left-hand side of equation (47) is just 2 ri c n (f). Using 
equations (51) and (52), equation (47) is rewritten as 


W > • 


~ 2t iff 

A< " ~ y + *» ' «, 


| exp + »)(/- tj] 


2*ic„(t ) = 



exp M« w + «)*J VO 


- exp ( - Kbt n - iT/2 - »/») (/ - fo)/rt)l| (59) 


+ Z Z v ,m AO 

^ ten m±n 


(53) 


The complex energy ? n for the intermediate state | n> is 
now defined as 


s e n + - ,T /Z 


(60) 


where 


Equation (47) then becomes 


vo * r — 

/-oo( w + 


exp [M/ 0 - 01 

+ + T-sH 


and 


(54) 


V' > V = 


K., exp (iw m /) 


X 


{ 


1 - exp (/((, - ?„) (f - 



/(/) 


.a j«f exp|i(u) + w nm +»a)r'| cji ) rf/'l exp ( - iu>r) < 

i»K/ + + ^r *2 ~ 


Equation (54) has two simple poles at 


(55) 


W 1 = — ( w m + '«) 


(56) 


and 


1 


£ £ ^/^mAO 

ten m*n 


(61) 


Equation (61) is the probability amplitude for the in- 
termediate resonating state. Note that /(f) comes entirely 
from the second term on the right-hand side of equa- 
tion (36). Thus, the first term on the right-hand side of 
equation (61) is the WICA contribution, and the second 
term is the contribution from higher order processes. 

As a check of the preceding results (eq. (61)). set /(f) 
= 0 and consider the matching of initial conditions in the 
WICA. With /(f) = 0, equation (53) becomes 




(57) 


VO = exp W“«I + in), ol VO 


(62) 
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which, for /< t ' yields 


c „<' < O = 0 


since / D (/</ 0 ) = 0 from equation (58). From equa- 
tion (61), for t-t 0 . 


c /.('o) = exp Mw„, + «)g / 0 (g = 0 (64) 
2*7# 


This expression satisfies equation (33) when equation (37) 
is substituted (after letting /— g. Thus, the initial condi- 
tions are satisfied in the WICA. 

The final-state amplitude is evaluated here by 
rewriting equation (5) as follows: 


= -U(') ex P ^ 
at in 


+ T I <)<(') exp (iw^r) (69) 


W = J° 


(w - (*>i) (-0) + w 2 ) 


using standard contour integration techniques. Thus, 
equations (63) and (64) satisfy the initial condition 
specified in equation (43). For t > t 0 and /(/) = 0, equa- 
tion (61) becomes 


where c n (i) is given by equation (61). Our interest lies in 
processes which proceed via the formation and decay of 
intermediate states. Thus, we will ignore the term involv- 
ing V kl , since it will be negligible or zero. Therefore, 
substituting equation (61) into equation (69) and in- 
tegrating from t 0 to t yields 

c kfO = £ ex P (“Vo> 




y„. exp (loV) 


x |l - exp li(«, - („) (t - f„)/#l| 


exp [/«*,(* - g] - 1 


exp l/(f* - 7 n ) ( t - g/#J - l 


The derivative of equation (66) is 


dc n« > to) V n, «P (*»„,') 


X ~ O ~ “*) 

x exp [i(t, - 7„) (t - g/#lj 


EE I Wm ' 


n l*n m±n 


f exp (iu kn t") I(t") dt" 


Equation (70) shows that the WICA is equivalent to 
neglecting at least third-order terms compared with 
second-order terms in the final-state amplitude. 

As was done previously for the intermediate-state 
amplitude, the final-state amplitude within the WICA is 
evaluated for the initial conditions, which are 


which, yields in the limit as t— t n 


c k (tst a ) = 0 


dc nU 0 ) ^WPf'Vo) 
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In equation (71), it is assumed that, before turning on the 
interaction at t Q> the system is in its initial state. Equa- 
tion (72) is obtained from equation (5), with the initial 
condition on the intermediate state specified by equa- 
tion (43) inserted. For /= f 0 , equation (70) satisfies equa- 
tion (71) trivially. Differentiating equation (70) and set- 
ting f = t 0 then satisfies equation (72). 

Finally, the integral /(/) is evaluated. Solution of 
equation (55) requires knowledge of c m (/'). Using only 
the lowest-order contribution 


If the result for a = 0 in equation (17) is compared with 
equation (75), the second-order and higher terms are 
identical if the real energies in equation (17) are replaced 
by complex energies ~ n . These results, taken with the 
similar results of references 9 and 10, demonstrate that it 
is permissible to replace real energies with complex 
energies, of the form given by equation (60), when study- 
ing resonances. 

Evaluation of Abrasion-Ablation T-Matrix 




(73) 


yields 


i' 


• a exp [i(u + m +ia)/'] tffvexp(-iaiO do) 


I r'“ 

/(')=i - — 
J-<* J/(u> + W 


- «)(-« [- ») 


(74) 

Evaluating equation (74) using the same contour integra- 
tion techniques used previously to solve l 0 (t) enable f he 
final-state amplitude from equation (70) to be written as 


c *(') = £ ~~ - ex P ('W,'< 


n*i e » ( i 
exp [i(t k - 


A 

t n ) (/ - g/A] - l) 

k ~~n ) 


exp [/w*,<r-g ]-l 


V kn V nl V !< e x P 

n l*rt t i 

( exp ( iu kl (t - g] - 1 

|( f * - «/) («/ - (/) (f, - f„) 
exp [io} k/ (t - gj - 1 

( f * - O (*/ - f /> («/ ~ 


exp [i(t k - T„) (/ - t Q )/fi] - 1 i 

- ~ n ) (e< ~ "«) < f / - ” n ) ) 

+ . . . (75) 


The basic Feynman diagram for projectile fragn u- 
tion, with fireball formation, is shown in figure 2. *he 
diagram is similar to those presented elsewhere (rels. 16 
through 18). Since our major interest is in the area of 
projectile fragmentation, target fragmentation and the 
formation and de-excitation of the fireball are not 
discussed in this paper. 

Recently, there has been considerable interest in 
relativistic coulomb dissociation and the excitation of 
nuclear giant resonances (refs. 19 through 25). The 
typical method for describing these excitations is through 
the interaction of the projectile with an equivalent target 
phonon. The excitation process is then described in a 
manner analogous to a photonuclear reaction. The basic 
Feynman diagram for coulomb dissociation is given in 
reference 20. The excitation of a giant resonance may be 
important in determining projectile prefragment charge 
dispersions (ref. 12) and requires further study. Since our 
primary interest is in projectile fragmentation, this pro- 
cess, rather than fireball formation, is emphasized in the 
Feynman diagram. If the fragmentation process is 
thought of as a collection of A projectile nucleons sud- 
denly being excited into a prefragment, then by analogy 
with coulomb dissociafion, the excitation process can be 
treated as the interaction between a phonon field and the 
initial A nucleons in the projectile. The Feynman 
diagram for this phonon excitation process is shown in 
figure 3. Since the fireball is difficult to depict, it is sim- 
ply represented as two separate projectile and target 
pieces. Figure 3 is only an alternvive, and more conve- 
nient, way of depicting the interaction shown in figure 2, 
and the diagrams yield identical results. In the actual 
fragmentation process shown in igure 2, there is no 
direct interaction between the target and excited 
prefragment. 

To simplify the phase space, the separate phase spaces 
of the individual particles can be replaced by a single 
phase-space factor describing the center of mass of those 
particles. Therefore, the fireball and target fragment are 
replaced here by target recoil T and projectile recoil R 
pieces, since we are not interested in the details of the 
phase spaces of these pieces. This Feynman diagram, 
which is used in this analysis, is shown in figure 4. The 
projectile recoil piece R is not the projectile prefragment 
P\ Figure 4 is exactly analogous to a 4-body final-state 
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Feynman diagram (Fig. 5) used to describe pion produc- 
tion via isobar formation and decay (refs. 9, 16, and 26 
through 30). This Feynman diagram has recently been 
replaced (refs. 16 and 31) by the 3-body diagram (fig. 6), 
in which the pion escapes but the nucleon is recaptured 
by the projectile nucleus. 

When developing expressions for the T-matrices, at- 
tention is focused on the interactions occurring in the 
projectile prefragment. (See Fig. 7.) As noted in the 
following section, the prefragment interactions in Fig- 
ures 4 and 7 represent only the lowest-order interaction. 
Higher-order terms are subsequently discussed in detail. 

The Two-Potential Problem 

In the basic fragmentation process it is assumed that 
the projectile nucleus experiences an interaction K 1 and 
undergoes abrasion. (See fig. 7.) The projectile prefrag- 
ment then experiences a different potential V 2 and 
ablates to yield the final fragment. This type of “two- 
potential” problem is considered in references 32 through 
34, and the simplest solution is stated in the pion- 
production work of Townsend et al. (ref. 30). 

The T-matrix expansion for abrasion-ablation (AA) is 

1% - V ki + VknG n ,V ni 


+ VknG„,V nm G mi V mi + . . . (76) 

where an Einstein summation convention on triply 
repeated indices is implied. The full interaction potential 
is separated into abrasion F 1 and ablation V 2 pieces as 

V = K 1 + V 2 (77) 


Inserting equation (77) into equation (76) and expanding 
yields 


T ** = v\. + V 2 + V 1 . G V X 

Kl ki ki kn ni ni 


+ V 2 . G V + k! G V 2 + V* G V 2 . 

kn m m kn m m kn m m 


+ v\ G y l G y . + v*. G y g v 2 . 

kn ni nm mi mi kn nt nm mi mi 


+ y[ g v 2 g k 1 + y[ g y 2 g v 2 . 

kn ni nm mi mi kn m nm mi mi 


+ v* g V g y' + v 2 . g y' g v 2 . 

kn m nm mi mi kn ni nm mi mi 


Each of these terms represents the higher-order 
generalizations of the Feynman diagram in figure 7. The 
Feynman diagrams, for the terms up to third order in 
equation (78), are shown in Figures 8 through 11, where 
the abrasion F 1 is represented by a wavy line (phonon) 
and the ablation V 2 is represented by a heavy solid line. 
As discussed subsequently, the Feynman diagrams permit 
the abrasion-ablation process to proceed in a time- 
reversed fashion with the ablation occurring prior to the 
abrasion. Also, these diagrams do not include exchange 
terms (ref. 35), because such exchange diagrams yield 
negligible contributions. 

Although equation (78) and its accompanying Feyn- 
man diagrams are quite complicated, they are shown only 
to third order and thus are not yet accurate enough for 
the present application. Ideally, interactions to infinite 
order should be included. This was effectively done for 
the abrasion step through the use of an eikonal scattering 
amplitude derived from an equivalent 1-body 
Schrodinger equation (refs. 2 and 4). This scatterng 
amplitude is simply and directly related to the abrasion 
T-matrix which, by definition, includes abrasion interac- 
tions to infinite order. Since the T-matrix in equation (78) 
includes both abrasion and ablation, the full abrasion 
T-matrix must be extracted to incorporate the previously 
developed abrasion formalism (ref. 4). In other words, 
the abrasion potentials, to infinite order, must be fac- 
tored out in equation (78). 

The technique for factoring out the full abrasion 
T-matrix is most clearly seen by again considering the 
Feynman diagrams. In figure 11, the abrasion process to 
all orders is depicted along with their corresponds 
mathematical expressions. The simplest approac> 
abrasion-ablation model would be to incorpora' 
first-order ablation into the infinite-order abra^on for- 
malism. (See Fig. 12.) This is done by simply adding abla- 
tion bubbles to the abrasion diagrams in figure 11. 
Mathematically, this is equivalent to simply multiplying 
the abrasion series by to yield 


V 2 G y 1 + V 2 G V 1 O K*. 

kn ni ni kn m nm mi mi 


+ v\ G V ,G, k! G y . + 

kn ni nl it Im mi mt 

= V 2 G (v l + K 1 G r 1 

kn ni \ ni nm mi mt 

+ v\gV, g y ] 

n.l h Im mi mi / 

= V 2 G T* br 
kn nt m 


(79) 


+ y g y 2 g V . + y g y 2 g v 2 

kn ni nm mi mi kn n: nm mi mt 

+ . . . (78) 


so that the factorization of the abrasion T-matrix is ap- 
parent. In figure 13, second-order ablation coupled with 
infinite-order abrasion is presented. The corresponding 
series is then of the form 
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+ (time-reversal terms) (80) 

where the time-reversal terms correspond to ablation oc- 
curring prior to abrasion. If these and any subsequent 
time-reversal terms from higher-order diagrams are ig- 
nored, the generalization to infinite-order abrasion- 
ablation clearly yields 



n 


This remarkable and extremely useful result arises entirely 
from equation (77) and the neglect of time-reversal 
Feynman diagrams. The neglect of these diagrams is 
henceforth dubbed as the “time-ordering approxima- 
tion. M Pilk^hn (eq. (3.15) of ref. 34) obtains a similar 
result when considering the time ordering of the time 
evolution operators (eq. (3.7) of ref. 34). For the purpose 
of this paper, equation (81) makes possible the incor- 
poration of the infinite-order ablation processes into the 
previously developed infinite-orde- abrasion formalism 
(refs. 2 and 4). 

Phase Space 

In this section the definitions and recurrence relations 
for the Lorentz-invariant, noninvariant, and “normal” 
density-of-states factors rre developed. 

Lorentz-invariant phase space is given by the 
restricted phase-space element (refs. 34 and 36) 


d Nips (e;p lf p 2 > . . P N ) 




(84) 


where v is the normalization volume. p or example, the 
corresponding recurrence relation for a 3-body phase 
space is given as 

d Nips (e;p lt p 2 ,p 3 ) 

= d Nips (e;/? c ,P 3 ) 


x d Nips (CftPit/^) dt c (85) 

Performing the momentum integrals over d Nips 
yields the usual density-of-states definition (ref. 37) 


f», (<*) m fd Nips (t k ; p x ,p 2 , . . p N ) 


= [jl-T'-' // 

L(2**)*J jj ■ 


L, 


x d i p ] d i p 2 




x d i p ] d*p 2 . . d*p N _ , (86) 


dLips (t\p x ,p 2 p N ) 


= (2ir)4 6 4 



1 " d% 

(2tt) 3N / a=l 2e i 


(82) 


where the 4-body recurrence relation, in terms of 2-body 
phase spaces, is 


where t k is the sum of the final energies of the N par- 
ticles. The recurrence relation for tfr N-body density of 
states is (ref. 16) 


. . n ) -/ * ’ * / ^ 2( f 12 * ' N ) P 2^\2 ‘ 


X ^2^ 12.. . N- 2) ■ ■ • ^2l f 12) ^ 12 df J2V &\2 . .yv~i 


where, for example, 


(87) 


d Lips U;p v p 2t p v p A ) = : d Lips (t>p c >p d ) 

(2tt)2 

X dUps(e d ;p\ y p 2 ) 

X dUp$(t c ;p 3 ,p 4 )dt c de d (83) 

where de c and dt d represent integrations over the “non- 
observed” particles. 

The noninvariant phase space, which differs slightly 
from that given on page 388 of reference 33, is defined as 


6 12 . . . N * e l + € 2 + • • • + *N (88) 

Abrasion-Ablation Model Derivation 

In previous work (refs. 6 and 38), abrasion-ablation 
cross sections have been determined by calculating abra- 
sion cross section: (refs. 3 and 4) which are then 
multiplied by an ablation probability obtained fre .1 com- 
pound nucleus decay probabilities (refs. 6 and 38). It is 
demonstrated in this section that this method of deter- 
mining abrasion-ablation cross sections arises solely f r om 
particular approximations to the general formalism 
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developed herein, and is therefore only a special case of 
this more general formalism. 

In terms of the transition rate, the total cross section 
is written as 

0 =-w (89) 

V 

where v is the incident velocity of the projectile. Inserting 
equations (30) and (86) into equation (89), the total 
abrasion-ablation cross section for the phase space 
associated with figure 4 is 

«*z) = hL -±L_ — 

fix (2 xff) 9 dt zxr T 

X /JT| T **\ ld ' p x d ' p r d ' p z (90) 

Using a recurrence relation (ref. 16) derived from equa- 
tion (86) 

Pa^zxrt^ = SSS 

x P2^Z\p) P Z)d d *ZX &ZXR <^I) 

demonstrates that <Pp x can be replaced by <Pp P in equa- 
tion (90) where 

<Pp p =<Pp zx (92) 


reasonableness of the classical argument. Equation (94) is 
sometimes justified quantum mechanically, especially 
when dealing with angular-momentum matrix elements 
(refs. 32 and 39) where theorems on Clebsch-Gordan 
coefficients are available (ref. 16, chapter 4). This is 
especially true, for example, for a single (one-level) reso- 
nant state involving several different angular-momentum 
projections M (ref. 39), where the summation over n 
simply becomes a summation over M for the single 
resonance state of a particular energy. This was also the 
case for the pion production work of references 26 
through 31, where there was only the single intermediate 
isobar £ resonance at a fixed energy but with various spin 
and isotopic spin projections. Norbury (ref. 16) has 
shown that equation *94) results from the spin-isospin 
Clebsch-Gordan algebra. Another example is the 
photonudear excitation of a compound nucleus where 
the formation of a resonant state of a single energy, but 
with different spin projections ‘ef. 40), justifies the use 
of the Bohr assumption when calculating (>,/?) cross sec- 
tions via compound nucleus formation and decay. In 
general, however, the preceding simplifications which 
justify the classical probability assumption do not hold 
for the abrasion-ablation process. For example, a par- 
ticular final projectile fragment could result from the 
ablation of numerous different prefragments, each with a 
quite different excitation energy. 

The partial width, which is simply a transition rate 
multiplied by Planck’s constant, is 


This, together with equation (81), allows the cross section 
in equation (90) to be written as 


o{Z) 


2tp v* d rrr y-abr 
Av (2xft) 9 <k ZX R T J” \n knm m 


X d*p p d*p T d*p z 


(93) 


= 2t 


v d 

(2^)5 *7 


/ 


n b J I 2 di Pz 


(95) 


Substituting equations (94) and (95) into equation (93) 
yields 


A major approximation is now introduced as 



*<Z) =£ 

n 


v v* d 

(2 xtf) 6 dt Z xRT' 


SSS rlG 


2 


x dtpp* d*p T dtp 


(96) 


which will henceforth be referred to as the “classical 
probability approximation,” since it involves the classical 
addition of probabilities (right-hand side) rather than the 
quantum mechanical addition of amplitudes (left-hand 
side). In essence, it involves ignoring the interference 
terms on the left-hand side of equation (94). It is our 
belief that the famous Bohr assumption for compound 
nucleus deca; (ref. 14), which justifies the separation of a 
two-step cross section (such as compound nucleus forma- 
tion and decay, or abrasion-ablation) into a product of 
formation and decay (partial width) cross sections, is 
based upon this classical probability approximation. The 
Bohr assumption is so widely used because of the 


which can be rewritten as 

”«>-££// rJcJICP 

n 

x dNips {ip RT \ p p „ p R .,p r )dt p , (97) 
The abrasion cross section is 

x d Nips (<p ‘ rt '<Pp ”Pr>Pt ) (98) 
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which yields 


If the branching ratio is defined as 


0{Z) = JF 2 * / r„ | G nj j ^ o n (A ) dtp. (99) 

n 

Inserting the Green's function, the abrasion-ablation 
cross section is 

<*z) = — £ f t ; a * (A) *p <>«» 

2r J +(r/2)2 

where the total T and partial widths are related by 

r=£r„ do!) 


To evaluate the integral in equation (100), the zero- 
width approximation (ref. 34) 


r-o (<„ - 1 ,) 2 +<r/ 2 p 


= *(«„ - *,) 


is introduced. Writing the energies explicitly as 


*n - e />' + c r + 

with an initial-state energy given by 


f * - 6 /> + e r 


and the Final-state energy as 


*k “ f X + f Z + *T + e H 


then conservation of energy 


fti-iyr (109) 

and is recognized as the usual ablation probability factor 
(refs. 6 and 38), then 

°( Z )='Z.Sn°nW 0 « 0 ) 


which is the standard abrasion-ablation cross-section 
result (refs. 6, 7, 38, 41, and 42). 

This result (eq. (’ 10)) can also be obtained from equa- 
tion (100) by an alternative method- Since o n (A ) is ob- 
tained by integrating over alt impact parameters, it is in- 
dependent of t p . . Taking it outside the integral enables 
equation (100) to be written as 

<■»> 

Inserting 


inside the integral in equation (111) and substituting 
equation (109) yields 

W* —■ *- (ll3) 

If g„ is independent of e p , (which merely requires T,, and T 
to possess the same energy dependence), then H can be 
taken outside the integral to yield 




” f / = + € z) 

Inserting equation (107) into equation (100) indicates a 
variable, intermediate, virtual resonance energy e p > 
centered about e x + e z , which is integrated over. The 
nature of the delta function in equation (102), however, 
destroys this quantum mechanical feature of virtual 
energy in the integral. The zero-width approximation, 
then, can be considered as another classical approxima- 
tion. Inserting equations (102) and (107) in equation (100) 
yields 

°(Z) = £(r„/rK(/i) (io«) 

n 


In principle, if the dependence of Ton t p . is known, then 
the integral can be calculated numerically if not 
analytically. If the zero- width approximation is inserted 
from equation (102), equation (110) is again obtained. 

Equation (1 10) is one of the central results of the pres- 
ent work. It represents a first-principles derivation of the 
usual abrasion-ablation cross section and results directly 
from: (1) the time-ordering approximation, (2) the 
classical probability approximation, and (3) the zero- 
width approximation. Clearly then, the most obvious im- 
provements to the abrasion-ablation theory would be to 
remove these assumptions (the time-ordering approxima- 
tion is the least important). The zero-width approxima- 
tion is removed in the next section. 
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Lorentz-invariant Differential Cross Section 

For an arbitrary number N of particles in the final 
state, the total cross section is 

. = ... / 

tlv L(2»a) , J dt k JJ J N-\ 

x\T A *\ 2 tPp\<flp 2 ...<Pp N _ l (115) 

For a specific final particle or fragment Z the differential 
cross section is 


^Z) _ 2TT f v XV i 

dP^ L< 2t, *) 3 J 


x 4 ~ /• - / / • • • / |r 4i4 | : 

dt k Z~\ Z+i ^V-ll ki I 


x d 3 p x . . d 3 p z _ , d 3 p z+ , . . . d 3 p N . , 


(116) 


This is formally cast into Lorentz-invariant form as 


where (ifio/dpP/t) is the Lorentz-invariant differential 
cross section obtained from equation (1 17). The notation 
“frame” in equations (119) and (120) denotes that the 
specified quantity is evaluated in the desired reference 
frame, which is completely arbitrary. The angular 
( do/dsi ) and spectral ( da/dt ) distributions are then ob- 
tained from equation (120) by performing the ap- 
propriate energy or angle integrations in that particular 
specified frame. 

For the abrasion-ablation Feynman diagram of fig- 
ure 4, the Lorentz-invariant differential cross section is 


d 3 a{Z) t z 2rr v3 d 

c i (dp i /t) z c 3 fly (2 xfl) 9 <k ZXRr 

xj/flEnxcr-'v “'or 

( 121 ) 

Employing the classical probability approximation and 
introducing dNips, equation (121) becomes 

<fia(Z) t z 2xv v 

c i {dp i /t) z ~ c 3 fly a*fi) 3 


dMZ) 

c 3 (dp 3 /f) z 

( Z [ 

"c 3 [ 

' V I rf J.. 

(2xA) 3 J dt k 



x/ .. 
Z+l 


■ d 3 p z , 


x 1 

• • • d 3 p s _ , 

(117) 


The which appears in equation (1 17) is cancelled by an 
e z which appears when the interaction matrix elements 
are explicitly evaluated. 

Using 

d*p = p 2 dpdn (118) 

it follows that 


*r/ra 2 K,rK?f 

n 

x dU\ps{( p RT ;p P ,p R ,p r ) (122) 

which, upon inserting equation (98), becomes 

d 3 o(Z) i z v 
c\dp 3 /t) z c 3 (2-irfi)’ 

X El^'l 2 \ G n,\ 2 °nW 023) 

n 

Inserting the Green's function yields 


d 2 o\ _ /p 2 \ 

dpd H / frame y f f frame dp 3 / 6 


(119) 


d 3 a{Z) ( Z v 
c\dp 3 / ( ) 7 ~ c 3 (2*f>) 3 


and 


= /*) -£±- 

fife d SlJ frame \^ 2 / f 


020 ) 


x 


E 


y-abl 2 
kn 




(e„ - €,) 2 + ( f/2) 2 


(124) 
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which demonstrates the need to evaluate the ablation 
(decay) matrix elements directly (when dea*ng with dif- 
ferential cross sections) rather than using partial widths 
(see eq. (110)), as was done with to'al cross sections. 

Equation (124) is a major result of the present work, 
since it provides the framework for calculating angular, 
spectral, and Lorentz-invariant distributions for the 
abrasion-ablation process in terms of the previously 
developed abrasion cross-section formalism (ref. 4). 
Fragment angular and spectral distributions can be ob- 
tained from equations (124) and (120) by performing the 
appropriate integrals. Clearly, the comparison of these 


distributions with experiment is of major importance in 
determining the accuracy of the present abrasion- 
ablation theory. Although these comparisons have not 
yet been done, equation (124) provides the framework 
necessary to do them. The remaining major task is to 
evaluate the ablation matrix element 7* bl . 

kn 

Langley Research Center 
National Aeronautics and Space Administration 
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Symbols 

A mass number of final projectile prefragment 

AA abrasion-ablation process 

c speed of light, 3x10* m/sec 

c,{r) initial-state probability amplitude 

c*(r) final-state probability amplitude 

c/^t) arbitrary-state probability amplitude 

c H (t) intermediate-state probability amplitude 

d Lips Lorentz-invariant phase-space element 

d Nips Noninvariant phase-space element, MeV " 1 
<Po differential cross-section element, mb 

dt energy element, MeV 

<Pp momentum phase-space element, (MeV/c) 3 

dll element of solid angle, sr 

F' highly excited fireball 

/„(«*)) Fourier transform of probability amplitude 

G Green’s operator, MeV - 1 

G m Green’s function, MeV~ 1 

g n branching ratio 

H full Hamiltonian, MeV 

H 0 unperturbed Hamiltonian, MeV 

ti Planck’s constant (6.58 x 10 -22 MeV-sec) 

7(0 time-dependent integral 

1 1 > initial-state vector 

| k> final-state vector 

M angular momentum projection, MeV-fm-c-* 

N number of bodies in final state 


|n> intermediate-state vector 
P projectile 

P(t) total transition probability 

P^t) single-state transition probability 

P' excited prefragment 

p momentum, MeV-c 1 

R recoil projectile 

T target 

T transition operator, MeV 

T k „ transition matrix element, MeV 

V recoil target 

/ time, sec 

t Q initial time, sec 

V interaction potential, MeV 

V kn matrix element of interaction potential, MeV 

V x abrasion potential, MeV 

V 2 ablation potential, MeV 

v incident projectile velocity, m/sec 

w transition rate, sec -1 

X ablated projectile particle* 

Z final projectile fragment 

z number of abraded protons 

a convergence factor 

T total decay width, MeV 

partial decay width, MeV 
7 photon 

A intermediate isobar resonance 
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energy shift, MeV 

Subscripts: 


Kronecker delta 

1 

initial state 


Dirac delta function on frequency, sec 

k 

final state 

€ 

energy, MeV 

i 

arbitrary state 

( i 

initial-state energy, MeV 

m 

arbitrary state 

( k 

final-state energy, MeV 

n 

intermediate state 

'* 

intermediate-state energy, MeV 

P 

projectile 

( z 

final energy of ablated particles, MeV 

T 

target 

m. 

( n 

intermediate complex energy, MeV 

Superscripts: 

V 

volume normalization element, fm J 

At 4 

abrasion-ablation 

*(<*> 

density of final states, MeV - 1 

abl 

ablation 

a 

cross section, mb 

abr 

abrasion 

* 

full wave function 



n 

solid angle, sr 



u 

angular frequency, sec - 1 



“kn 

frequency difference between states | k> and 
|n>, sec -1 
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Figure 3. Feynman diagram of equivalent target phonon excitation of projectile 
prefragment. 



20 


Figure 4. Equivalent Feynman diagram (lowest order; no time 
reversal) of projectile prefragmentation used in this work. 





Figure 7. Lowest-order projectile prefrag- 
ment interactions. Wavy line represents 
abrasion V 1 , and heavy solid line 
represents ablation V 2 . 




Figure 8. First-order matrix elements of abrasion-ablation T-matrix 
in equation (78). 
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Figure 9. Second-order matrix elements of abrasion-ablation T-matrix in equation (78). 
(The first abrasion process on left-hand bubbles is that of figure 4. All other abrasion 
processes are only indicated schematically. Particle number conservation is not shown 
explicitly as it is for ablation.) 
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V, 1 G .V 1 G .V 1 . V, 1 G V 1 G .V 2 . V, 1 G G .V 1 . 

kn m nm mi mi kn ni nm mi mi kn n mi mi 



V 1 G V 2 G V 2 V 2 G V 1 G V 1 V 2 G V 1 G V 2 

kn ni nm mi mi kn ni nm mi mi kn ni nm mi mi 



V 2 G V 2 G V 1 V 2 G V 2 G V 2 

kn ni nm mi mi kn ni nm mi mi 


Figure 10. Third-order matrix elements of abrasion-ablation T-matrix in equation (78). (Abrasions 
are indicated schematically as in fig. 9.) 
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V 1 G V 1 G V 1 
nl li lm mi mi 


Figure 1 1 . Abrasion to all orders. 



V 2 G V 1 G V 1 G V 1 
kn ni nl li lm mi mi 


Figure 12. Abrasion to all orders with ablation to first order. 
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Figure 13. Abrasion to all orders with Ablation to second order. 



